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Abstract 

For two nearby pulsars, formed close to the galactic plane, their host 
open stellar clusters have been found. 



1 Birth site for pulsar B 2224 + 65 



Cordes et al. (1993) discovered the complex consisting of pulsar B 2224 
+ 65 and of the gaseous nebula "Guitar", the pulsar being located at 
the end of its "neck" . The pulsar is specified by a great proper motion 
which may be due to its proximity to the Sun or, otherwise, it is a very 
peculiar object even among the pulsars. The only procedure used so 
far to measure the distance using the dispersion measure of the pulsar, 
DM = 35.3 ps/cm~ 3 , yields D ~ 1.95 kps. There is, however, a 
ring of gas around the pulsar ionized by a shock wave (the density 
jump is characterized by a factor of 10 12 (see Kaplan & Pikelner, 
1979)). Electrons which contribute to the dispersion measure of the 
pulsar may partly belong to this gaseous ring. Therefore, the value 
Do ~ 1.95 kps seems to be overestimated. The discovery of soft X- 
ray radiation from the nebula "Guitar" , that contains the moving pulsar 
(Romani et al., 1997), also provides evidence for the complex to be close 
to the Sun. A more precise distance to the pulsar could be estimated 
from its trigonometric parallax when available. 

Characteristic pulsar age is t c = 0.5 P/P = 1.12 Myr (Taylor et 
al., 1993) where P is the pulse period, P, its first derivative in seconds. 
Characteristic ages of pulsars are found to be systematically high. A 
real age is determined more precisely using the expression for so-called 
"kinematic age" (Smith, 1977): 

t k = 0.5r D ln(l + * c /(0.5td)). 

The parameter of magnetic field attenuation is estimated by us (Pskovskii 
& Dorofeyev, 2001) to be td = 1-8 ± 0.2 Myr. Thus, kinematic age of 
pulsar 2224+65 is t k = 0.74 ± 0.10 Myr. 

The position angle PA of the pulsar motion has been preliminarily 
found to be 52° which indicates that the pulsar has arrived at its current 
site on the celestial sphere from the Illrd quadrant, that is, from the area 
within a(19 h + 22 h ), 5(+70° -=- 65°). To determine the site of its birth 
more accurately one needs to take into account systematic variations 
in pulsar radial velocity and proper motion, though its spatial velocity 
may be regarded as being constant. 

If the pulsar has been moving towards the Sun for the past 0.7 
Myr, it is likely to have started from the region close to an open stellar 
cluster localized not far from the Sun in the same quadrant. A possible 
candidate is cluster M39 (NGC7092). 
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2 Formulae for calculation of the pulsar 
starting site 

A trajectory of motion for a pulsar having a constant spatial velocity is 
rectilinear, but its projection on the celestial sphere is an arc of a great 
circle. We shall use the rectilinear equatorial reference frame XYZ with 
the origin O at the point of observations (Fig. 1). The current position 
II(a , S , 1950) of the pulsar belongs to the plane YZ. It corresponds 
to the position U(x , y , z Q ) in the rectilinear frame, where x Q = 0, 
y = D cos S , z = D sin J . The components of the pulsar spatial 
velocity V are directed, correspondingly 



V a (V a , y , z ) — parallel to the axis OX, 



V$(0, Do — Vs smS Q , z + Vs cosS Q ) — along the meridian, 
V r (0, y + V r cos S , z + V r sin S ) — along the ray II O. 

These components are obtained in km/s if expressed in terms of ob- 
servational values V a = 4.74£> 0) u a , V 5 = 4.74 D /j, s , V t = 4.74 D 0/ u (the 
tangential pulsar velocity for the current epoch), V r = V cos(\ Here \x, 
H a , n s are the proper motion and its components along the axes OX 
and OY expressed in arc seconds per year and corrected for the motion 
of the Galaxy and for the motion with respect to the Sun; ( is the an- 
gle between the rotational axis of the pulsar and the direction towards 
an observer. The expression V r = V cos( follows from the result es- 
tablished earlier (Tademaru, 1977; Pskovskii & Dorofeyev, 1987); the 
directions for both V and Vt are supposed to be known from the "rule 
of signs" for ( (Pskovskii & Dorofeyev, 2001). 

The system of equations for a straight line parallel to the vector 

V(V a , V r cos S — V§ sin S Q , V r sin S Q + Vs cos S ) 

is given by 

x y -y Q z - z 



V a V r cos S — V5 sin S a V r sin S + Vs cos S Q 



0.978t fc = tf, (1) 



where x, y, z are the rectilinear coordinates of a pulsar at a starting 
moment (in parsecs), 0.978 = fifth is a dimensional factor in units of 
ps/km/s. By their sense, the values tk and fi are negative. As the 
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value D is not known, and neither are the velocity and its components 
expressed in terms of D 0} the system ([I]) can be transformed as follows: 

x'/Do = 4.74tf/z a , 

y'/D = (l±4.74^cotC)cos5 -4.74^sin^, (2) 

z'/D = (1 ±4.74$// cot sin S Q + 4.74^ cos S § , 

{D'/D f = {x'/D f + {y'/D f + {z'/D ) 2 . 

D 1 is the pulsar distance at a starting moment. The equatorial coordi- 
nates for the starting moment a', 5' are given by: 



, y' ■ ( x '\ 

a = ao H tj arcsm — 

cos o \D' J 



z' 



5' = arctanf— ). 

y 

The coordinates x' M ,y' M ,x' M ;a' M ,8' M for cluster M39 at the moment 
— tk are obtained using the formulae ([[]). 



3 The birth site of pulsar B 2224+65 

The parameters of pulsar B 2224+65 and of the open stellar cluster 
M39 which are currently observed, as well as those calculated for the 
moment of birth, are shown in the second and the third columns of 
Table 1, respectively. The other columns of the Table contain findings 
for the objects to be reviewed below. The nomenclature of the rows is 
explained in the text. The data for the pulsars are given according to 
Taylor et al. (1993), and those for the open stellar clusters, according 
to Rastorguev & Glushkova (1999). The angle $ has been determined 
by means of our approach (Pskovskii & Dorofeyev, 1998). 

Since the distance D Q to pulsar B 2224+65 is not known, the values 
of x' I D Q , y 'I D Q , etc. are to be calculated. For the rest of the pulsars and 
for all the stellar clusters the values x', y', etc. have been determined 
as the distances D a are available for them. 

Because of significant errors involved in the technique of determining 
initial coordinates of a pulsar using its current estimated values fi, V r , 
D Q , and one can anticipate, in the best case, but the full (or a partial) 
superposition of an error box of a pulsar and that of a stellar cluster. 
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In our case, the distance to pulsar B 2224+65 estimated at the moment 

of birth cannot be matched with that for the stellar cluster M 39. 
Table 1 



Objects 
cuo (h, m) 

son 

C'/yr) 
C'/yr) 
n C'/yr) 

*(°) 

K) 
D (pc) 

Vt (km/s) 

V r (km/s) 

V (km/s) 

t k (Myr) 

t c (Myr) 

•d 

(x'/Dq) , x' (pc) 

(y'/Do) , y' ( P c) 

(z'/Do) , z' (pc) 
(D'/D ) , D' (pc) 

5>n 

Aa' (h, m) 
a' (h, m) 



PSR B2224+65 
108.6, +6.8 

22 h 2A m 

+65.34 
0.146+0.003 
0.113+0.003 
0.185+0.003 

41+3.0 

(200 ± 50) 
(200 ± 30) 
(-230 + 30) 
(-300 ± 40) 

1.12 
0.74+0.010 
+0.74+0.010 
(-0.49 + 0.07) 
(1.13 -=-0.07) 
(1.41 + 0.010) 
(1.87 + 0.09) 

51+11 
l h 07 m ± 32 m 



M39 


PSR B192929+10 


NG( 


92.5, -2.3 


47.4, -3.9 


36.1 


21 /l 30.4 m 


19 /l 29.9 m 


18 h 


+48 .22 


+ 10.88 


+ 


=0.0504+0.0013 


0.096+0.006 


0.000 


-0.016+0.003 


0.050+0.004 


0.002 


0.0191+0.0013 


0.108+0.006 


0.002 




24.4+3.0 





1.5 






350+100 


170+40 


38C 


32+7 


90+20 


4: 


-7.5+5 


-190+50 


-2 


33+7 


210+50 


3( 




3.09 






1.34+0.16 






-90+30 




12+6 


400+80 


0: 


220+90 


30+50 


410 


280+110 


410+80 


40 


360+100 


3.7+20 


410 


53+33 


3.7+2.0 


1.4 


+ 12 m ±9 m 


54 m ± 18 m 


+ 


21 /i 45 m ± gm 


lS^™ ± 18 m 


18^25.5 



21^17 m ±32 m 

Therefore, the knowledge of possible uncertainties in a' and 8' is of 
great importance. To estimate them, the value of D' /D Q was suggested 
to possess a relative error which corresponds to the "bb" class in the 
classification proposed by Taylor et al. (1993), that is, 46 per cent. 
Errors of the values in the formulae have been suggested to be mean- 
square ones deduced from expressions for corresponding functions. 

The error box for the stellar cluster M 39 has been calculated taking 
into account the angular diameter of the cluster having an extended 
corona with a diameter of ~ 1.5° (Artiukhina, 1970; Barkhatova & 
Pylskaya, 1978). 

Fig. 2 shows the error boxes of these objects to be partly superim- 
posed at the starting moment of the pulsar. The pulsar therefore can be 
considered to be at the distance D' at the moment in question, this dis- 
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tance corresponding to cluster M 39. This fact allows us to calculate the 
values D Q , V, Vt, and V r for the pulsar. The findings are shown in Ta- 
ble 1 in brackets. The current distance to the pulsar, D Q = 200 ± 50 ps, 
is obtained to be dramatically small, whereas kinematic parameters are 
shown to have rather common values. A dispersion measure of the pul- 
sar in such a case is mainly due to electrons in a shock wave rather than 
to interstellar medium. 

A spatial trace of the open stellar cluster M39 (Fig. 2) is far shorter 
than that of pulsar B 2224+65. The distance to M39 has remained 
practically the same, but the directions of pulsar and cluster motions 
are almost opposite. 

Taking into account that M39 is presumably a member of a kine- 
matic group in Ursa Major (Lloyd Evans & Meadows, 1964; Eggen, 
1965), an error box of the pulsar starting position proves to lie com- 
pletely inside a field of the group. Unfortunately, it is impossible to 
determine kinematics of the object before the start, hence we cannot 
establish with certainty whether it belongs to the group or not. 



4 The birth site of nearby pulsarl929+10 

The coincidence of the birth site of pulsar B 2224+65 with the field of 
the stellar cluster M 39, however, does not exclude the possibility that 
some single pulsars may be produced outside stellar clusters. A num- 
ber of young single pulsars, namely B 0531+21, B 0833-45, B 0611+22, 
B 0656+14, etc., which have not travelled very far from their starting 
points are not related to any stellar clusters. Some pulsars are found 
to start at high galactic ^-coordinates (Pskovskii & Dorofeyev, 2001). 
From statistical findings correlations were shown to exist among young 
pulsars and stellar associations (Mdzinarishvili, 1997). In addition, pul- 
sars discovered in globular stellar clusters belong to a quite specific 
group, as well as millisecond pulsars and those in binary systems. 

Let us consider other pulsars which are close to us at the present 
time and which have well established estimates of proper motion, the 
distances D , and ages within the limits 0.7 < % < 1.5 Myr. In ad- 
dition to pulsar B 2224+65 considered above, these conditions are met 
for the two more pulsars, B 1929+10 and B 1133+16, the latter having 
unprecedented proper motion. Pulsar B 1133+16 is believed to have 
been born at a high galactic latitude, both for V r > and for V r < 0, 
so it has no host stellar cluster. 
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Unlike pulsar B 2224+65, determinations of the birth sites of the two 
pulsars mentioned above are based on distances established rather well 
by Taylor et al. (1993). The positions of these pulsars at a starting 
moment have been calculated using formulae (1). 

The data concerning a birth site of the pulsar and the parameters 
of its host cluster NGC 6638 are given in the corresponding columns of 
Table 1. The value of ( has been estimated by our technique (Pskovskii 
& Dorofeyev, 1998). Cluster NGC 6633 has low proper motion, its error 
box lying almost entirely within that of pulsar B 1929+10 at the moment 
of its start (Fig. 3). 



5 Conclusion 

It has been shown that two nearby pulsars, which started from the 
galactic plane, may have been related, at a moment of their birth, to 
open stellar clusters. It is worthwhile pointing out that the both pulsars 
were reported to be weak sources of X-rays (Wang et al., 1993; Cordes 
et al, 1993). 

There is good reason to believe that actually there are more pul- 
sars related to open clusters. This can be deduced from the work of 
Mdzinarishvili (1997) containing evidence for the correlation of young 
pulsar (t c < 2 Myr) distribution close to the galactic plane with that of 
O-associations in the vicinity of the Sun. 

The authors wish to express gratitude to the researches of the As- 
tronomical Department of the MSU Physics Faculty, Professor A. S. 
Rastorguev and Assistant Professor Ye. V. Glushkova, who have given 
the data for the open stellar clusters. 
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Figure captions 



Fig.l. Projections of a pulsar radial velocity V r and of the compo- 
nents V a and V$ of a pulsar tangential velocity Vt on axes in the rectilin- 
ear equatorial reference frame X, Y, Z. The case of V r > is assumed. 
The following designations are used: Ayi = cos Jo, A^i = V^sin^o, 
A 2/2 = V r cos5o, Az2 = V r smSo. 

Fig. 2. Positions of pulsar B 2224+65 (empty squares) and of the 
open stellar cluster M 39 (empty circles) in the sky at the present time. 
For the pulsar starting moment the objects are indicated correspond- 
ingly by filled squares and filled circles which lie in the centres of their 
error boxes. The dashed line stands for the error box of the cluster; 
the continuous line, for the pulsar. The heavy line signifies the galactic 
equator. Coordinates are given for the epoch 1950.0. 

Fig. 3. Positions of pulsar B 1929+10 and of the open stellar cluster 
NGC6633 in the sky. The designations are the same as to Fig. 2. 
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